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We report a detailed analysis of the Drude weights for both thermal and spin transport in one 
dimensional spin- 1/2 systems by means of exact diagonalization and analytic approaches at finite 
temperatures. Transport properties are studied first for the integrable XXZ model and second for 
various nonintegrable systems such as the dimerized chain, the frustrated chain, and the spin ladder. 
We compare our results obtained by exact diagonalization and mean-field theory with those of the 
Bethe ansatz, bosonization and other numerical studies in the case of the anisotropic Heisenberg 
model both in the gapless and gapped regime. In particular, we find indications that the Drude 
weight for spin transport is finite in the thermodynamic limit for the isotropic chain. For the 
nonintegrable models, a finite-size analysis of the numerical data for the Drude weights is presented 
covering the entire parameter space of the dimerized and frustrated chain. We also discuss which 
conclusions can be drawn from bosonization regarding the question of whether the Drude weights 
are finite or not. One of our main results is that the Drude weights vanish in the thermodynamic 
limit for nonintegrable models. 



I. INTRODUCTION 

Due to their relevance for modeling the magnetic 
properties of several quasi-one-dimensional materials, 
dimerized and frustrated spin-1/2 chains and ladders 
are models of great and current interest. Recently, 
growing experimental evidence has been found that 
magnetic excitations can contribute significantly to 
the thermal conductivity of various quasi-one and 
quasi-two-dimensional materials 1-11 . Stimulated by 
these observations, many theoretical activities have 
addressed the issue of heat transport in one-dimensional 
spin systems 12-22 . While spin transport has been a topic 
of numerous theoretical investigations 22-45 , the theory 
of thermal transport is less well understood. 
One of the key questions is to understand under which 
conditions transport is ballistic, i.e., dissipationlcss. 
The criterion for this is the existence of a singularity 
at zero frequency in the real part of the conductivity. 
Therefore, one is interested in the integrated weight of 
this singularity - the so called Drude weight - in the 
thermodynamic limit. The appearance of a nonzero 
Drude weight is often ascribed to the influence of 
conservation laws on transport properties 30 ' 38,43 ' 46 ' 47 . 
For example, in the case of the Heisenberg chain the 
energy-current operator is conserved 30,48 implying a 
nonzero thermal Drude weight at all temperatures. An- 
other widely discussed and related issue is the difference 
between transport in integrable models compared to 
nonintegrable ones 15,17 ' 20 ' 27 ' 33 ' 41,43_ 

The purpose of the present paper is to provide a system- 
atic study of both the Drude weight for spin and thermal 
transport at finite temperatures by means of exact diag- 
onalization on finite systems and analytic methods. The 
dependence on exchange coupling anisotropy, frustration 
and dimcrization will be clarified. 

We start with an overview of the results of previous 



and related works. For the anisotropic Heisenberg chain, 
it is possible in principle to compute thermodynamic 
quantities exactly with the Bethe ansatz at arbitrary 
temperatures. The Drude weight for thermal transport 
has been obtained in the antifcrromagnetic regime of 
this model along this line 14,49 . Numerical studies 15 ' 17 
have provided the thermal Drude weight in the antifcr- 
romagnetic regime, in the gapped, ferromagnetic phase 
and for the isotropic, ferromagnetic chain. Here, we 
extend on such analysis by adding data for the gapless, 
ferromagnetic regime. 

In a first numerical work devoted to the issue of 
thermal transport in nonintegrable models, Alvarez and 
Gros 15 have conjectured that the Drude weight is gener- 
ically finite in dimerized and frustrated spin systems 
although the energy-current operator is not conserved 
in these cases. However, we have argued in Refs. 17,21 
that this conclusion cannot be sustained for gapped, 
frustrated chains if larger system sizes are included in 
the finite-size analysis. In this paper, we extend our 
parameter study of the thermal Drude weight to include 
dimerized chains and spin ladders also. The main result 
is that the numerical data are best interpreted in terms 
of a vanishing thermal Drude weight in nonintegrable 
systems. 

Recently, however, the thermal Drude weight has also 
been computed by means of analytic approaches yielding 
a finite Drude weight at low temperatures 18 ' 20 . In Ref. 
18, the spin ladder and the dimerized XY model have 
been studied by mapping to non-interacting models. 
While this is exact in the latter case, it is an approxi- 
mation in the case of the spin ladder. For instance, the 
influence of incommensurate umklapp-scattering terms 
on transport properties of massive models is not yet fully 
understood. Bosonization was applied in Ref. 20 to the 
cases of the dimerized and the frustrated chain leading 
to the interesting result that certain umklapp-scattering 
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terms do not spoil the conservation of the energy current. 

Spin transport in spin-1/2 models is equivalent 
to charge transport of (spinlcss) fcrmions and an enor- 
mous amount of work has been devoted to this field 22-45 . 
The situation of spin transport in the integrable model, 
i.e., the Heisenberg chain, is different from the case of 
thermal transport since the spin-current operator is 
only conserved in the case of free fermions (XY model). 
Nevertheless, the Drude weight is finite in the gapless 
regime while numerical 29,32 and analytical 36 studies 
have found indications that it may vanish in the gapped 
cases. The case of the isotropic chain is still the subject 
of discussion 31,33 ' 37 - 41 ' 43 ' 44 . 

At present, no final agreement about the results 
from different Bcthe ansatz computations 34,37 ' 50 ' 51 for 
the Drude weight has been achieved. For example, the 
Bethe ansatz computations by Zotos 37 predict the Drude 
weight of the isotropic chain to be zero while Kliimper 
and co-workers 50 have found a finite Drude weight in 
this case. The latter conclusion is in agreement with 
some numerical works - exact diagonalization 31 ' 33 (ED) 
and Quantum-Monte Carlo simulations 41 (QMC) - as 
well as analytic approaches 43 . Very recently, however, 
Long et al. 44 have developed a novel Lanczos method 
for finite temperatures which has been applied to spin 
transport in the anisotropic Heisenberg chain with up 
to N = 28 sites. They interpret their numerical results 
in terms of a vanishing Drude weight for the isotropic 
chain. 

Regarding the temperature dependence of the Drude 
weight in the gapless regime contradicting results are 
reported in the literature 37,41 ' 43 ' 50 ' 51 . 

As far as nonintegrable models are concerned, the 
effect of frustration on spin transport has been studied 
at zero temperature 26 , while at finite temperatures, 
Ising-like interactions with distance up to 3 have been 
investigated with ED 29 ' 33 and QMC 41 . Fujimoto and 
Kawakami have discussed transport in nonintegrable 
models in the low-energy limit with bosonization. They 
find a nonzero Drude weight provided that particle-hole 
symmetry is broken in agreement with Ref. 30. The 
physical reason is the finite overlap of the spin-current 
operator, which is not conserved in the presence of any 
umklapp scattering, to the energy-current operator. 
In their work on charge transport in a one-dimensional 
system 38 , Rosch and Andrei have argued that the Drude 
weight is zero in the presence of dangerously irrelevant 
perturbations since then all relevant conservation laws 
are broken. On the basis of this result, a zero Drude 
weight is naturally expected in all nonintegrable and 
gapless spin models. 

It will be an additional focus of the present paper to 
provide a complete as possible analysis of spin transport 
in nonintegrable models such as the dimcrized and 
frustrated chains as well as spin ladders. As a main 
result, the numerical data indicate a vanishing of the 
Drude weight in the thermodynamic limit consistent 



with Refs. 29,38. 

Although it is beyond the scope of this paper to 
give any conclusive explanations of recent transport 
experiments on low-dimensional spin materials, we will 
provide a list of materials for which this work could 
be of relevance. Whenever possible we will make con- 
nection to these experiments and discuss implications 
of our results for their interpretation. First, there 
are surprisingly large magnetic thermal conductivities 
observed in the compounds (Sr,La,Ca) 14 Cu2404i 2 ~ 5 . 
Here, a strong magnetic contribution to the thermal 
conductivity is believed to be mediated by magnetic 
excitations of spin ladders. Second, similar findings 
have been reported for the spin chain materials SrCu02, 
Sr2Cu03 6 and BaCu2Si20y 7 . Third, we mention the 
anorganic Spin-Peierls material CuGcOa 1 although the 
physical nature of the heat carrying excitations is still 
under controversial discussion 1 ' 52 . It has been suggested 
that the magnetic properties of this material can be 
described to some extent by a dimcrized and frustrated 
chain. Finally, strong evidence for thermal transport 
mediated by magnon excitations in two-dimensional 
cupratc compounds has also been found 8-11 . 

The plan of the paper is the following. In Sec. II, 
we introduce the model which is a dimerized and frus- 
trated Heisenberg chain including an exchange-coupling 
anisotropy. Next, we summarize the basic relations 
for the transport coefficients and, particularly, for the 
Drude weights as they follow from linear response theory 
in Sec. III. The definitions for the current operators 
will be given at the end of this section. Our results are 
presented in Sees. IV and V. First, we discuss transport 
properties of the anisotropic Heisenberg chain using 
mean-field theory and exact diagonalization. The results 
will be compared to Bethe ansatz computations 14 ' 37,49 , 
QMC 41 , and other numerical studies 15,32,33,44 as well 
as analytical works 39,43 . Second, we consider several 
nonintegrable models in Sec. V. Using bosonization, 
we argue that vanishing Drude weights are generically 
expected in massless 38 and massive nonintegrable spin 
models. Numerical data will be shown for both thermal 
and spin transport and will be related to the results of 
other groups 15,18,20,43 . Finally, our main conclusions are 
summarized in Sec. VI. 



II. MODEL 

In this work we consider frustrated and dimerized spin- 
1/2 chains as sketched in Fig. 1. The Hamiltonian reads: 

N 

H = ^(A,a,A) (1) 
i=i 

hi(A,a,X) = J[Xihi,i + i(A) + ah l ,i+2{ A )} ( 2 ) 

A M+i (A) = (#^+H.c.)/2 + ASfS? +i . (3) 
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FIG. 1: Sketch of a frustrated and dimerized 5 = 1/2 chain. 
The limiting cases are: a = 0, A = 1: the XX Z chain; (ii) 
a > 0, A = 1: the frustrated chain; (iii) a = 0, < A < 1: the 
dimerized chain; (iv) a > 0, A = 0: the two-leg spin ladder. 



hi(A,a,X) from Eq. (2) is the local energy density and 
N the number of sites. We set A; = 1 if I is even and 
A; = A otherwise. In Eq. (3), i = 1 for nearest- and i = 2 
for next-nearest-neighbor interactions. Throughout this 
work, periodic boundary conditions are used. All quanti- 
ties are measured in units of the exchange coupling con- 
stant J where J is positive. Mainly, we focus on the 
following limiting cases of this model: (i) a = 0, A = 1, 
the XX Z chain; (ii) a > 0, A = 1, the frustrated chain; 
(iii) a = 0,0<A<l, the dimerized chain; and (iv) 
a > 0, A = 0, the spin ladder. 

We discuss various values for the anisotropy in the inte- 
grable case (i.e., a = 0, A = 1) including the XY model 
(A = 0). In the case of the two-leg spin ladder (A = 0), 
a = Jii / J± is equal to the ratio of the leg coupling Jn to 
the rung coupling J±. 

It is instructive to write down the Hamiltonian in terms of 
fermionic operators. This will later be used as the start- 
ing point for the mean-field theory in Sec. IV A. By per- 
forming a Jordan- Wigner transformation (see, e.g., Ref. 



53), the spin-1/2 operators S l 
spinless fermions 



±,z 



can be mapped onto 



g+ =e <T*« c t 



(4) 



„(t) 



destroys(creates) a fcrmion on site I. The string op- 
erator <I>; reads $; = Yli=i c l c i- For the Hamiltonian 
from Eq. (1), we arrive at (see e.g., Ref. 54) 



III. KUBO FORMULAE FOR « AND a 

In this section we briefly review the underlying equa- 
tions for the thermal and the spin conductivity and 
we define the energy-current operator jth and the spin- 
current operator j s . 

A. Drude weight for heat and spin transport 

The Kubo formulae for the thermal conductivity k(uj) 
and the spin conductivity cr(u) read 53,55-57 



oo 



-i{ui+iQ + )t 



dT(j t h[s](-t - ir)j th [ s ]}. 

(6) 

The expressions in brackets [■] refer to spin transport; 
the index 'th' denotes thermal transport and 's' spin 
transport. (3 = 1/T is the inverse temperature while 
the brackets ( . ) denote the thermodynamic expectation 
value. The real part of k(u>) and cr(u>), respectively, can 
be decomposed into a ^-function at lu = and a regular 
part 

ReK[a](cj) = D th[s] (T)S(u) + Re K[a} Tes (cj) (7) 

Here, we are interested in the weight D t h[s] (T) of the 5- 
type singularity at u) = 0. A nonzero value of D t h[s](7 1 ) 
corresponds to a divergent conductivity whereas a van- 
ishing _D th [ s ](T) implies a normal conducting behavior 
if there is a contribution from the regular part of the 
conductivity at low frequencies or insulating behavior if 
ReK[cr] rcg (w = 0) = (see Ref. 25). 

In the remainder of the paper, the dependence on sys- 
tem size of the Drude weight will be explicitly denoted 
by D t M s ](N,T) while D t M s ](T) refers to the thermody- 
namic limit N — * oo. 

A spectral representation of /t[cr](cj) in terms of eigen- 
states H\n) = E n \n) yields an expression for both Drude 
weights 30 (Z is the partition function) 



dL [s] (n 7 t) 



ZN 



E 



|(m|j th[s] |n)| 2 . (8) 



H = 



[jM(c;Vi+H.c.) 
i L z 

+A(c]c l cj +1 c l+1 -4c l + -)} 
+ a {\i4 c i+2 +H.c.)(l -4 +1 c l+1 ) 
+A(cjqcj +2 q +2 - 4 c l + j)} 



(5) 



Note the presence of correlated hopping terms, i.e., 
c i c i+2 c i+i c i+i +H.c, in Eq. (5). Such contributions arc 
absent in the usual standard extended Hubbard models. 



Note that the exponent £ t h[i 
chosen according to 



in the prefactor needs to be 



tth = 2 or t s = l, (9) 

i.e., we expect D t h ~ 1/T 2 and D s ~ 1/T at high tem- 
peratures. 

For spin transport, we will numerically analyze Eq. (8) 
and a second expression which can be derived from linear 
response theory 23 ' 27 



D s "(iV,T) 



7T 

N 



-f) 



2 



E 



-0E, 



E„ — E,, 



(10) 
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The operator f is related to the kinetic energy. Note 
that in the derivation of Eq. (10), the /-sum rule (see 
Ref. 23 and references therein) is exploited. It is now 
instructive to make a connection to Kohn's formula for 
the Drude weight D s to illustrate the relation between 
Dl(N,T) and £> s 7/ (TV, T) . As a by-product, we will ob- 
tain the appropriate definition of T. 
Equation (10) is equivalent to the generalization of 
Kohn's formula 58 for finite temperatures 27 



D[\T) = 



ZN 



-pE n 



d 2 E n (<i>) 



a$2 



(ii) 



$=0 



To arrive at this expression, one considers finite rings of 
length TV being pierced by the flux <f>, i.e., a static twist 
about the z axis. The flux-dependent Hamiltonian H(&) 
with exchange couplings J r ^ being nonzero for r < r 
reads 26 

N r 1 

= EEM (^ W ^+Hx.) 

1=1 r=l 

+AS?S? +r }. (12) 
Expanding if ($) up to second order in <f> yields 23,26 



if ($) = if ($ = 0) 



TV 



3s 



27V 2 



T . 



(13) 



The expression for the current operator j s will be given 
in Sec. IIIB where we derive it from the equation of con- 
tinuity. In terms of spin operators, T is given by 



N ro 



r=EE^ r (^\+H,c, 



(14) 



1=1 r=l 



For the models investigated in our paper, we have ro < 2. 
In the fermionic picture (see Eq. (5)), T is the sum of 
the contributions to the kinetic energy from both the 
nearest-neighbor and next-nearest-neighbor hopping and 
also the correlated hopping terms c[, jC^cj, 2 c ( + H.c. 
The appearance of the extra factor of 4 in front of the 
term proportional to J2 = a J, which arises from the 
factor r 2 in Eq. (14), can be understood by considering 
the limit of two decoupled chains {J a = const; J — ► 0) 
of length TV/2. Rescaling the size-dependent prefactors 
in Eq. (13) on N/2 cancels the factor of 4 and results in 
twice the Hamiltonian of a single chain with N/2 sites. 
By applying second-order perturbation theory to H(&) 
it is straightforward to arrive at Eq. (10) (see, e.g., Ref. 
23). 

The relation between Df (TV, T) and £> 7/ (TV, T) has been 
pointed out in Refs. 25,30,33,35. While 

D S (T)= lim D I S {N,T)= Urn L> b 7/ (TV,T) 



values on finite systems at T = (sec Rcfs. 24,35 for 
examples) whereas £> 7 (TV, T) > is strictly fulfilled by 
evaluation of Eq. (8) (see, e.g., Ref. 30 for details). Note 
in this context that the difference D**(N, T) - D 7 (TV, T) 
is proportional to the second derivative of the free energy 
with respect to the flux $ in the limit of vanishing mo- 
mentum. In general, this quantity measures the super- 
fluid density 25 ' 30,35 in more than one spatial dimension. 
However, for a one-dimensional model, it vanishes in the 
thermodynamic limit 30,35 . 

Eq. (8) can also be related to Kohn's formula. While the 
latter one is the sum of the curvature of energy levels, 
Dg(N, T) equals the sum of squares of first derivatives of 
the energy levels with respect to the flux <i> 33 



Dl{N,T) 



ZN 



£< 



.-(IE, 



$=0 



(15) 



In this paper, Eqs. (8) and (10) will be used to com- 
pute the Drude weights while Eqs. (11) and (15) have 
been given for completeness. We will show in Sec. IV B 3 
that the difference between D\{N ,T) and D^(N,T) is 
negligibly small at sufficiently high temperatures already 
on finite systems. At low temperatures, only D^ 1 (N,T) 
yields a good description of the temperature dependence. 

If the current operator j t h[s] is a conserved quantity, 
Eq. (8) can be rewritten as 

DL [S] (N,T) = Z^L^e-^^.jIn) (16) 



ZN 

n 

N \Jth[i 



(17) 



implying that a static expectation value has to be com- 
puted which is a considerable simplification. Analo- 
gously, we obtain 



Dl\N,T) = -{-f). 



(18) 



from Eq. (10) under the condition [H,j s ] = 0. The 
main goal of this work is to investigate whether T?th[s] 
has a finite value in the thermodynamic limit (TV — > 00). 
Typically, the finite-size dependence is well controlled at 
high temperatures T> T where all states \n) are occu- 
pied with equal probability e~@ n — > 1. In this case, we 
see from Eq. (8) that 



lim \T tth ^ 

T^oo 



dLJn,t)] 



C th[s] (N) (19) 



N 



N- 



with a prefactor C t hf! 



Cth [s] (TV) 



Z^N 



IHith[s]|")| 5 



(20) 



holds for all temperatures T > in the thermodynamic 
limit, the two expressions are not equivalent on finite 
systems. For instance, Eq. (10) may result in negative 



In the limit T — > 00, the partition function Z 
equal to the dimension of the Hilbert space. 
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B. The current operators 

The current operators obey the following equations of 
continuity 

d t hi = i[H, hi] = -(jtv+i -jth.O (21) 
d t S? = i[H,Sf]=-(j s , l+ i-j St i) (22) 

where hi is the local energy density, Sf the local magneti- 
zation density and Jth[s],i are the local current operators. 
From Eqs. (21) and (22), we obtain solutions for the cur- 
rent operators for arbitrarily long-ranged interactions. If 
[hi± r , hi] 7^ holds only for r < r , the energy-current 
operator reads 17 ' 48,59 

TV TV r a - 1 

jth = £ jth, I = *£ £ [hl-r-i, hi +n ]. (23) 

1=1 1=1 n,r=0 

The spin current is given by 23 ' 26 

TV TV i-o-l 

Js = £ j S ,i = *£ £ fo-r-l. S l+nl (24) 
1=1 1=1 n,r=0 

with [hi_ r ,Sf] for r < tq. In this paper, we re- 
strict ourselves to tq < 2. The same result can be ob- 
tained from Eq. (13). Independent of the model, the 
local magnetization density is given by Sf whereas we 
use Eqs. (2) and (3) as a definition for the local energy 
density. For the XY model, both quantities jth and j s 
commute with the Hamiltonian, i.e., [H , jth[s]] = 0. 
While [H xxz , jth] = remains valid for arbitrary A 
and a = 0; A = 1 30 > 48 ; the conservation of j s is immedi- 
ately broken once the anisotropy A has a nonzero value. 
Dimerization and frustration which spoil the integrabil- 
ity of the XXZ model also lead to [H, j t h[ s ]] (except 
for A = a = 0). 

IV. XXZ MODEL 

In this section we discuss transport properties of the 
XXZ model. Since this model is solvable via the 
Bethe ansatz, certain thermodynamic quantities like 
susceptibility 60,61 or specific heat 14,61 are accessible with 
this technique. The thermal Drude weight D t h has re- 
cently been computed along this line by Kliimper and 
Sakai 14,49 for arbitrary A > 0. For spin transport, finite- 
temperature Bethe ansatz computations have been per- 
formed by both Zotos 37 and Kliimper and coworkers 50 . 
However, these two groups find different results for the 
temperature dependence which is ascribed to conceptual 
problems of computing the Drude weight from a single 
microstate at a given temperature 51,62 . 

First, we apply mean-field theory based on the Jordan- 
Wigner transformation and compute the Drude weights 
in the gapless, antiferromagnetic regime < A < 1. This 
approach yields a good approximation to the exact re- 
sults from the Bethe ansatz (compare Rcf. 17) for the 



thermal conductivity whereas for spin transport, it fails 
to describe the A-dependence of the Drude weight D s . 
Second, we present our results from exact diagonaliza- 
tion focusing mainly on spin transport. 

A. Mean-field theory 

In this section, we consider an approximate treat- 
ment of the Hamiltonian of the anisotropic chain within 
mean- field theory. First, this approach will provide ex- 
act expressions for the Drude weights in the case of free 
fermions (A = 0) which is useful for a check of the numer- 
ical implementation. Second, wc will see that the thermal 
Drude weight of the XXZ chain can be well described 
by this approximation for |A| < l 17 which is, however, 
not the case for spin-transport. 

Within mean-field theory, the interaction terms in the 
fermionized Hamiltonian from Eq. (5) are simplified by 
means of a Hartree-Fock decomposition. Since details of 
the procedure can be found in the literature (see, e.g., 
Ref. 54), we prefer to quote directly the result for the 
mean-field Hamiltonian (k labels the momentum) 

H M F = E e fe4 c fe- ( 25 ) 

fc 

Hmf is diagonal in momentum space with the tight- 
binding dispersion 

e fe = -t(A,T)cos(k). (26) 

The hopping amplitude i(A,T) := -(1 + 2A(T)A) is 
modified via the one-particle expectation value 

A(T) = (cjc l+1 ) = - f dk cos(fc)/( efe ), (27) 

* JO 

which needs to be determined self-consistently. Note that 
/(e) = l/(exp(/3e) + 1) is the Fermi function. 
In the mean-field description, the current operators are 
simply given by 

j t r=$>w is MF = £^ nfc! (28) 
k A,- 

with velocity Vk = dek/dk and n% = c|.c fc . The spinon 
velocity, i.e., vff F = v k (T = 0), follows from Eq. (26) 
v™ F = 1 + 2A(T = 0) • A with A(T = 0) = 1/tt. From 
Eq. (17), we obtain the Drude weights 

D T(T) - ^£( e ^) 2 / 2 ( e <0 e£fc/T > (29) 
fc 

D™ F (T) = ^£ v fc7 2 ^K fc/T (30) 
fc 

In the case of free fermions, i.e., A = 0, Eqs. (29) and 
(30) are exact and can be evaluated both on finite sys- 
tems and in the thermodynamic limit. This provides a 
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useful check for the numerical implementation; in partic- 
ular since the spin-current operator j a is independent of 
A (sec Eq. (24)). 

For thermal transport, the mean-field approach pro- 
vides a good approximation to D t h(T) for |A| < 1 as we 
have shown in Ref. 17 and as can be seen in the inset 
of Fig. 2 (a) where we compare Dth(T) from the mean- 
field approach (dotted line) with the Bcthe ansatz results 
from Ref. 14 (solid line) for A = 0.5. 
Regarding the spin transport, we note that the mean- field 
theory does not give a qualitatively correct description of 
D B (A). Equation (30) results in D S MF (T = 0) = < IF (A). 
Basically, the Drude weight increases with A here be- 
cause D^ IF ~ (— T) essentially measures the increase of 
the bandwidth with A. This is in contrast to the exact 
expression for D S (T = 0) 23 

D 8 (T = 0) = ^-^L A = cos( 7 ) (31) 
4 7(7r — 7) 

where D s decreases with A due to the presence of the 
interaction terms in Eq. (5). The reason for the failure of 
the mean-field theory is that the current operator j s 
is conserved for < A < 1 which is not true for the 
full Hamiltonian (1) and the respective current operator 
(24). Note that the exact Drude weight takes the form 
D s (T = 0) ~ K vq where K is the Luttingcr parameter 
and vq the velocity at T = (see, e.g., Ref. 39). It is 
precisely the factor of K that is missing in the mean- 
field result. 



B. Numerical results for the XXZ model 

Now we turn to the discussion of the results from ex- 
act diagonalization for (i) the thermal Drude weight D t h 
and (ii) the Drude weight D s for spin transport of the 
anisotropic chain. Before going into details, let us out- 
line the structure of this section and point out the main 
results: 

First, we discuss technical aspects of the numerical pro- 
cedure which will also be of relevance for the nonintc- 
grable models in Sec. V. In particular, the importance of 
degenerate states and their numerical identification will 
be emphasized. 

Second, the numerical results for the thermal Drude 
weight D t h(N,T) for A = ±0.5,2 will be presented and 
compared to the Bethe ansatz 14,49 . In addition, the high- 
temperature prefactor Cth is computed and discussed as 
a function of A. 

Third, a detailed parameter study of the Drude weight 
D s is performed. Results for D^' U (N,T) are shown 
for A = 0.5, ±1,-2 as well as for A > 1. We dis- 
cuss the temperature dependence of both D^(N,T) and 
Dg 1 (N, T) and point out that they are indistinguish- 
able at high temperatures, but exhibit a completely dif- 
ferent finite-size scaling 33 at low T. This result, i.e., 
D*(N,T) = D S /7 (7V,T) for T large enough, is of impor- 
tance since the Drude weight could in principle also be 



obtained analytically from Eq. (8) or (15), respectively, 
which has to our knowledge not been attempted so far. 
Further, the influence of logarithmic finite-size correc- 
tions at T = for the isotropic chain is mentioned 39 . 
Finally, and most important, we analyze the finite-size 
scaling of the high-temperature prefactor C s and its de- 
pendence on A. While we can unambiguously con- 
firm that limjv->oo C s ( N) > for |A| < 1 in agree- 
ment with the results of other groups 29 : 32 ' 33 * 37 . 41 . 43 , 4 ^ 
the data for A = 1 clearly indicate a finite Drude weight 
D S (T > 0) as well. The latter observation is in agreement 
with Refs. 31,33 ' 41,43 ' 50 but contradicts the conclusions of 
Refs. 37 - 44 . 



1. Technical remarks on the numerical procedure 

We start with several technical remarks on the numer- 
ical procedure which arc relevant for both integrablc and 
nonintegrable models. We have performed complete di- 
agonalization for chains with N < 18 sites exploiting con- 
servation of the z-component S* ot = ^ z S? of the total 
spin, translational invariancc, and spin-inversion symme- 
try in the S z = subspaces of systems with even N. The 
latter symmetry is respected by the energy-current oper- 
ator jth but not by the spin-current operator j s . The 
dimensions of the largest subspaces for a given momen- 
tum k arc w 2400 for 5 t z ot = 1 and w 2700 for 5f ot = 
at TV = 18. In the latter case, the dimension is almost 
reduced by a factor of 2 by spin-inversion symmetry for 
the subspaces with odd and even sign under this symme- 
try. 

Another important aspect is the identification of de- 
generate states (i.e., E n = E m ) in subspaces labeled by 
S^ ot and momentum k. This is necessary both in the 
evaluation of Eqs. (8) and (10) but becomes irrelevant if 
the respective current operator is conserved, leading to 
the simpler expression, Eq. (16). The latter is possible 
for thermal transport in the XXZ model. 
For spin transport, however, we have [H,j s ] ^= for 
A / 0. The (integrated) distribution of level spacings 
AE n is shown in Fig. 3. There, the number /(e) of level 
spacings AE n = E n+ \ — E n ,E n < E n+ \ of adjacent en- 
ergy levels being smaller than a given value of e is plotted 
versus e for A = 0.5 

J ( e )= E E !■ ( 32 ) 

(S*, t ,fc)AE„<e 

It is sufficient to analyze all subspaces with given Sf ot 
and momentum k separately and sum over all subspaces 
thereafter (indicated by the first sum in Eq. (32)). 
The spectrum displays some characteristic features: first, 
the value of /(e) for e — > oo equals the dimension of the 
Hilbert space 2^ minus the number of subspaces (S* ot , k). 
Second, a large fraction of degenerate states is present 
and third, the integrated distribution of level spacings is 
constant for 10 -9 J < e < 10~ 6 J for the system sizes in- 
vestigated here. This suggests that adjacent energy levels 
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FIG. 2: Thermal Drude weight of the XXZ chain for: A = 
0.5 (panel (a)), A = —0.5 (panel (b)) and A = 2 (inset in 
panel (b)). Dashed(dot-dashed) lines: ED data for even(odd) 
N. Solid lines: Bethe ansatz (BA) results from Refs. 14, 
49. Thin, solid line in (a): high-temperature limit D t h(T) — 
Cth/T 2 from Ref. 14. Dotted line in the inset of (a): mean- 
field (MF) approximation. 

are typically separated by AE n < 10 -9 J if they are de- 
generate. 

This separation allows for an identification of degenerate 
states by imposing the following criterion in the numer- 
ical analysis: energy levels with \E n +i — E n \ < e cu t = 
1CP 8 J arc degenerate. By evaluation of Eq. (8) for the 
thermal Drudc weight we find agreement with Eq. (16) 
and the Bethe ansatz 14 for this choice of e cu t but signifi- 
cant deviations at temperatures T ~ J if larger or lower 
values for the cutoff energy are used. 

2. Thermal transport in the XXZ model 

Numerical results for the thermal Drude weight 
Dth(N,T) have previously been presented for A = 1 
in Refs. 15,17 and for A = -1; -2; 10 in Rcf. 17. Wc 
enhance on the latter study by adding numerical data 
for the gapless antifcrromagnctic and the ferromagnetic 
regime which is shown in Fig. 2 (panel (a): A = 0.5, (b) 
A = —0.5). Dashed (dot-dashed) lines display ED-data 
for N even(odd). Since we are using periodic boundary 
conditions, the twofold degeneracy of the ground state 
for chains with an odd number of sites leads to a diver- 
gence of D th (N,T) for T -> 0. The Drudc weight of 
systems with an even number of sites is characterized by 
an exponential suppression D^[N,T) ~ e~ AFS / T a ^ \ ow 
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FIG. 3: Distribution of level spacings in the spectrum of finite 
chains with 8 < N < 18 (bottom to top as indicated by the 
arrow) for A = 0.5 (solid lines: N even, dashed lines: N odd). 
AE n is the difference of adjacent energy levels in subspaces 
classified by total S* ot and momentum k. The number 7(e) 
of AE„ with AE n < e summed over all subspaces is plotted 
versus e (see Eq. (32)). 

temperatures where Aps is the finite-size gap of the sys- 
tem. 

A striking feature of D th (N, T) of the XXZ model is the 
rapid convergence with N at high temperatures. The dif- 
ference between the numerical data for the largest sys- 
tem investigated for A = 0.5 (i.e., N = 18) and the 
Bethe ansatz curve 14 (solid line in Fig. 2 (a)) is smaller 
than 10~ 2 J 2 for temperatures T > 0.23 J. Basically, the 
convergence of the numerical data is as good as in the 
case of other thermodynamic quantities such as specific 
heat Cy and susceptibility \- We stress that the relation 
Dth{T) ~ Cy(T) only holds at high and low tempera- 
tures (see Ref. 14). 

The antiferromagnetic, gapped regime has been discussed 
in Ref. 17. We mention that in the meantime Bethe 
ansatz computations 49 have been extended to A > 1. 
The numerical data are in agreement with these results 
(see inset of Fig. 2 (b): solid line , BA; dashed line, ED 
for N = 18 at A = 2). 

We continue by a discussion of the high-tempcraturc 
prefactor Cth defined in Eq. (19). As is evident from Fig. 
2 (a) and (b), the value of Cth := limjv_oo C t h(N) can 
be numerically determined already from systems with a 
comparably small number of sites (e.g., N = 12) with 
very good accuracy. Here we are interested in the A- 
dependence of C t h as shown in Fig. 4 (squares: ED). 
Note the excellent agreement with the analytic expres- 
sion from the Bethe ansatz 14 (solid line in Fig. 4) 

C^A) = ^(3 + S ^^m . (33) 
64 \ sm(arccos(A)) / 

For A = (i.e., free fermions), we get Cth (A = 0) = 
7T J 4 /32 which can also be obtained from Eq. (29). 
There are two main features of Cth (A): (i) Cth (A) = 
Ct h (-A) (sec Fig. 4 (a)) and (ii) lim A ^ 00 [C th (A)/A 4 ] = 
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FIG. 4: High-temperature prefactor Cth for the XXZ model 
in the gap less (left panel) (|A| < 1) and gapped, antiferro- 
magnetic regime (right panel) (A > 1). The numerical data 
(ED) are in perfect agreement with the exact expression for 
C t h from the Bethe ansatz 14 (BA). 



(see Fig. 4 (b)). In Fig. 4 (b), C t h is measured in units 
of A. 

The first property can be seen from the following obser- 
vations: changing the sign of A is an antiunitary trans- 
formation and essentially turns H(—A) = ~H(A) and 
E n (~A) = —E„(A) while the eigenvectors remain un- 
changed (this follows from an additional rotation by tt 
about the z axis on all sites with even site index). This 
transformation leaves the energy-current operator un- 
changed, i.e., jth(~ A) = jth(A). From Eq. (20), we can 
then conclude C t h(A) = C t h(— A). 

For J/A -> 0, the Ising limit H = J2i SfSf +1 is ap- 
proached where the only possible excitations are local 
spin flips. Here, no current can flow and, consistent with 
this notion, [hi 7 hi + i] = leading to jth = (see Eq. 
(23)) and D th = for J/A = 0. 

3. Spin transport in the XXZ model 

In the following, a survey of the results for D s for ar- 
bitrary values of the anisotropy will be given. First, the 
gapless, antiferromagnetic regime (0 < A < 1; Fig. 5) 
will be discussed with a particular focus on the isotropic 
chain (A = 1). Second, we comment on the finite- 
size data for D s for the gapped, antiferromagnetic case 
(A > 1; Fig. 6 (a), (b)), and third, results for the ferro- 
magnetic regime will be shown (A < 0; Fig. 6 (c), (d)). 
Finally, the dependence of the high-temperature prefac- 
tor C s on both anisotropy and system size is analyzed 
(Fig. 7). 

Numerical results for D^ ,H (N,T) in the gapless, an- 
tiferromagnetic regime are shown in Fig. 5 for A = 0.5 
(panel (a)) and A = 1 (panel (b)). Note that our re- 
sults for D ! s (N, T) agree with the data for A = 0.4 (not 
shown in the figures) and TV < 14 published in Ref. 33 
by Narozhny et al. 

First, we concentrate on the case of A = 0.5. In panel 
(a), we compare data from Eqs. (8) (D^(N,T), dashed 
lines) and (10) (D^ 1 (TV, T), solid lines) confirming that 
these two expressions are equivalent at high temperatures 




1 • 1 ■ 1 1 1 • 1 1 1 

0.2 0.4 0.6 0.8 1 

T[J] 

FIG. 5: Spin transport: Drude weight for A = 0.5 (panel (a)) 
and A = 1 (panel (b)). Dashed lines are numerical data ob- 
tained for D' a {N, T) from Eq. (8); solid lines denote D"(iV, T) 
from Eq. (10). Arrows indicate increasing system size. Note 
that Dl 1 (N, T) shows a better convergence with TV at low 
temperatures compared to Dl(N,T). The inset shows the 
Drude weight at T = for A = 0.5 (diamonds) and A = 1 
(squares) as a function of the inverse system size. The dashed 
lines are fits applied to the subsets with even (solid symbols) 
and odd (open symbols) N . Open circles at T = denote the 
exact values for D S (T = 0) 23 . 



(T > 0.5 J, depending on system size). 

At low temperatures, D{(N,T) shows much slower 
convergence with TV and essential features of the tem- 
perature dependence are only present in the data for 
Df / (iV,T) which arc the finite value at T = and the 
vanishing slope of D'J (TV, T) for T -> 0. The latter obser- 
vation (i.e., dDg I (N, T)/dT = at T = 0) is consistent 
with the Bethe ansatz by Klumper and co-workers 50,51 . 
In general, the functional form of D S (T) at low T for 
|A| < 1 is 

D S (T) =D s (T = 0)- const ■ T a (34) 

where the exponent a depends on the anisotropy. In Ref. 
37, the expression a = was derived (The integer v 
parameterizes the anisotropy via A = cos(tt/is)). Tak- 
ing this result, Eq. (34) would imply D S (T) = D S (T = 
0) — const ■ T for A = 0.5 which is obviously not con- 
sistent with our numerical data. For |A| < 1 and low 
T, the temperature dependence seems to be described 
by the expressions derived by Fujimoto and Kawakami 43 
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which arc compatible with QMC and the ED presented 
here. 

We have also determined D z 1 (N, T = 0) numerically 
by evaluating Eq. (10) in the subspace containing the 
ground state (see inset of Fig. 5). Using Eq. (11), one 
could go to larger systems than N = 18 since only the 
curvature of the ground state is needed. However, the 
main topic of this work is the Drude weight at finite tem- 
peratures while the finite-size corrections for the T = 0- 
Drude weight have been computed in Ref. 39. 
Data for £>f 7 (A, T = 0) are plotted versus 1/N in the in- 
set of Fig. 5 for A = 0.5 (diamonds) and A = 1 (squares). 
The data from finite systems with an even number of 
sites form a monotonically decreasing sequence with N at 
T = and small temperatures (see Fig. 5 (a), T < 0.2 J) 
while the data for odd N are a monotonically increas- 
ing sequence. Thus, the results for D^ I (N,T) and even 
./V provide an upper bound and those from systems with 
odd N a lower bound for D S (T) at low temperatures. 
From Ref. 39, the leading finite-size corrections at T = 
are available 

D II (N,T = 0)=D s (T = 0) + ^ I + .., (35) 

with [i = 2 for A < 0.5. Performing fits according 
to Eq. (35) at T = separately for even(odd) TV re- 
sults in D S (T = 0) = 0.9747(0.9717) J for A = 0.5 
which is in very good agreement with the exact result 23 
D S {T = 0) = 0.97428J. 

At the isotropic point, i.e., A = 1, the curves dis- 
play similar features as for A = 0.5: (i) a vanishing 
slope for T — > 0, (ii) a monotonic decrease at high 
temperatures, and (iii) D I S I (2N,T) < D S /7 (2A - 2,T) 
(D I S I (2N + 1, T) > D I S I (2N — 1, T)) at low temperatures 
T < 0.1J. 

However, the finite-size data at the isotropic point and 
T = seem to follow £> 7/ (A, T = 0) = A + B/N in con- 
trast to the case of A = 0.5 as can be seen in the inset 
of Fig. 5 (b). For A, we find A f=a 0.847J which com- 
pares well with numerical results obtained by the Lanc- 
zos method reported in Ref. 26. Admittedly, a good ap- 
proximation to the exact value of D S (T = 0) at A = 1 
for N — > oo cannot be obtained from the numerical data 
since the system sizes are far too small (see Fig. 5 (b)). In 
fact, from the work of Laflorencie et al. 39 it is known that 
the relevant and leading finite-size correction at T = 
and A = 1 is a logarithmic term. This is because umk- 
lapp scattering is a marginally irrelevant perturbation in 
this case. Similar to the susceptibility x( T ) 39 ' 60 , D a (T) 
is expected to show a sharp drop for T — > accompanied 
with a diverging slope at T = 50 in the thermodynamic 
limit. 

At sufficient large temperatures, we believe that the nu- 
merical data for A = 1 presented in Fig. (5) (b) give 
the correct picture of the temperature dependence of the 
Drude weight. However, Fujimoto and Kawakami 43 have 
recently obtained an analytic expression for D S (T) in 
the low-energy limit with conformal field theory which 
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FIG. 6: Spin transport. Panels (a) and (b): Drude weight 
D^(N,T) as a function of A. (N = 9, 13, 17; in (a): T = 1J; 
(b): T = 0. Panel (c): Drude weight D* Z (N,T) for A = -1 
and 8 < N < 18 (N even, top to bottom). Panel (d): Drude 
weight for A = —2 and = 12, 14, 17 (top to bottom; dashed 
lines: D£(N, T) Eq. (8), solid lines: D S "(7V,T) Eq. (10)). 
Arrows indicate increasing system size. 



is compatible with our numerical data for |A| < 1 but 
not for A = 1. Here, Fujimoto and Kawakami find 
D S (T) < D S (T = 0) while the data shown in Fig. 5 (b) 
seemingly support the opposite relation. Despite this dis- 
crepancy, our results do nevertheless support the notion 
of a finite D S (T > 0) at A = 1. This is substantiated 
by the analysis of the high-temperature prefactor C S (N) 
(see Eq. (20)) as we will discuss in detail below. 

For larger A (i.e., A > 1), the monotonic increase 
of £>g 1 (N, T) at low temperatures for odd N changes 
to a monotonic decrease. This is illustrated in Fig. 
6 (b) where D* Z (N,T = 0) is plotted versus A for 
N = 9,13,17 (compare Refs. 39,63). The cross-over in 
the monotony occurs at A m 1.2, i.e., in the gapped 
regime. Since D II (N,T) w const for small T, the be- 
havior at T = is characteristic for the low-temperature 
regime. At larger temperatures, D 11 (N,T) is a mono- 
tonically decreasing function for both even and odd N 
(see panel (a) of Fig. 6 for odd N and T = 1J). 

Regarding the ferromagnetic regime (i.e., A < 0), we 
concentrate on A = — 1 and A = — 2. The results for 
A = -1 plotted in Fig. 6 (c) indicate D TI (N,T) « const 
at low T with D S (T = 0) « 0.523(5)J. However, since 
the low-energy spectrum for A = —1 is of comparable 
complexity as for A = 1, one may expect nontrivial finite- 
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size corrections which could lead to a different tempera- 
ture dependence at low T. If the latter is true, then the 
system sizes are too small to draw conclusions about the 
behavior of D I S I {N, T) at very low T. 

The Drude weight in the gapped, ferromagnetic regime 
is expected to show a behavior analogous to that of 
A > 1. For instance, D^(N,T) is monotonically de- 
creasing with N at all temperatures irrespective of odd- 
even effects. Interestingly, (N, T) and D^(N, T) turn 
out to be indistinguishable for A < — 1 and N large 
enough which is illustrated in Fig. 6 (d) . This plot shows 
Dl' n (N,T) for N = 12, 14, 17 at A = -2 where dashed 
lines denote (TV, T) from Eq. (8) and solid lines stem 
from Eq. (10) (D^(N, T)). 

To conclude this section, we discuss the high- 
temperature prefactor C s both as a function of anisotropy 
and system size. In Fig. 7 (a), C S (N) is shown versus 
l/N for A = 0, 0.5, 0.6, 1, 1.5. Results for the case of free 
fcrmions (A = 0) have been deduced from Eq. (30) on 
finite systems in the fermionic picture. To fix the abso- 
lute values we note that C s — ir/8 in our notation for the 
free-fermion case A = 0. 

As is evident from Fig. 7 (a), C S (N) roughly follows 



Cs(JV) 



b 

N 



(36) 



The same finite-size extrapolation has been used for 
Dl (N, T) ■ T in Rcf. 33 for even-numbered systems and 
N < 14 at T = 50 J. However, by comparing the direct 
computation of C S (N) from Eq. (20) with fitting 



Dl' H (N,T) 



C S (N) Ci(JV) 



T 



T 2 



(37) 



to Dl' n (N, T), wc find that, often, more than one term 
in Eq. (37) needs to be taken into account to recover 
the result for C S (N) from Eq. (20). Consequently, it is 
preferable to compute C S (N) directly from Eq. (20). 
The data can be well extrapolated according to Eq. (36) 
for A = 0.5 and A > 1. A subtlety arises from the strong 
differences between even and odd-numbered systems for 
intermediate values of A such as A = 0.6 (see Fig. 7 (a)). 
In these cases (i.e., A = 0.2, 0.4, 0.6, 0.8, 0.9) we have sep- 
arately performed fits to the subsets with even and odd 
N. C s is then estimated by averaging the results from 
the fits. The large error bars for A = 0.2, 0.4, 0.6, 0.8, 0.9 
are due to these large differences in the extrapolated val- 
ues of subsets with even and odd N. 
Following this procedure, we obtain C s = (0.15±0.03) J 2 
for A = 0.6 while the extrapolation of the even-numbered 
systems yields C s = (0.123 ± 0.001)J 2 . The latter 
value compares well to the data published in Rcf. 33 
(C s = (0.119 ± 0.004) J 2 ). The strong finite-size effects 
between even and odd numbered systems, however, in- 
dicate that additional finite-size corrections apart from 
a simple 1 /iV-term must become relevant for larger N if 
the sequences of even and odd N converge to the same 
value in the thermodynamic limit. A difference in the 
extrapolated value for large N of these two subsets does 
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FIG. 7: Spin transport. Panel (a): high-temperature pref- 
actor C S (N) versus l/N for A = 0, 0.5, 0.6, 1, 1.5 (solid(open) 
symbols: even(odd) N). The dashed lines are fits according to 
Eq. (36) (see text for details). Panel (b): extrapolated high- 
temperature weight C a (A) as a function of the anisotropy A 
(diamonds) and the data for N = 17, 18 (squares, circles). 
The lines are guides to the eyes. 



not seem to be plausible. 

Note that it has been suggested in Refs. 32,44 to partly 
include spectral weight from low frequencies of Recr(<x>) 
to obtain C s for all those A which cannot be written as 
A = cos{tt /v),v integer. However, within such a pro- 
cedure, a criterion is necessary to decide up to which 
frequency one should integrate Reer(a;) and it goes obvi- 
ously beyond the usual definition of the Drude weight via 
Kohn's formula (11). In this paper, we prefer to confine 
our analysis to Eq. (10), or (11) respectively. 

At the isotropic point A = 1, we have fitted the nu- 
merical data for C S (N) using Eq. (36) while varying the 
range of data points (A m ; n < N < N max ). Stable fit re- 
sults are obtained for 11 < N < -/V max with N max > 14 



and for N„ 



< N < 18 with AT min < 14. Wc find 



C s = (0.020±0.006) J 2 indicating a finite Drude weight in 
the thermodynamic limit. This result is slightly smaller 
than the data reported in Rcf. 33 (C s = (0.046±0.005) J 2 
there for N = 6, 8, ... , 14) but the latter is regained if 
we use our data for ./V = 6, . . . , 14. Since we observe 
C S (A,N) = C S {-A,N), this implies D S (T > 0) > for 
the isotropic, ferromagnetic chain also. 

The results from the extrapolation are summarized in 
Fig. 7 (b) together with our data for N = 17 and N = 18. 
The plot suggests that (i) D S (T) > for A = 1 and (ii) 
D S (T) = for A > 1.5 within our numerical accuracy. 
Respective conclusions can be drawn for A < — 1 since 
C S (A) = C s (— A). Wc stress that for intermediate Ising- 
like anisotropics (i.e., 1 < A < 1.5) the system sizes may 
still be too small for an unambiguous confirmation of the 
conjecture 29 ' 32,36 D S (A > 1) = 0. In particular, the pos- 
sibility of a finite Drude weight in the gapped regime 
cannot be ruled out on the basis of the numerical data 
even though the Drude weight is zero at T = 0. An 
example for such a scenario (i.e., D S (T = 0) = but 
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D S (T > 0) > 0) has been given in Ref. 35. 
Very recently, Long et al. 44 have applied a newly de- 
veloped finite-temperature Lanczos method to compute 
C S (N) for N = 24,26,28. Compared to our data, their 
results strongly deviate from the fits to Eq. (36) for 
systems with N < 18 for all values of A presented in 
Fig. 7 (a). 



V. NONINTEGRABLE MODELS 

In this section we address the issue of transport in non- 
integrable models by means of bosonization and exact di- 
agonalization. As discussed in Sec. IIIB, both the spin- 
current and heat-current operators are not conserved in 
the presence of frustration and dimerization except for 
the case of a dimcrized XY chain 18 . 

The original conjecture by Zotos and co-workers 27-29 
stated that the Drude weights are expected to vanish in 
nonintegrablc models and we will argue in the following 
that our exact diagonalization study corroborates this 
statement. 

In a first numerical work by Alvarez and Gros 15 on ther- 
mal transport, the data obtained by complete and exact 
diagonalization of systems with N < 14 have been inter- 
preted in favor of a nonzero Drude weight D t h(T > 0) > 
in the cases of spin ladders and frustrated spin chains. 
However, we have argued 17,21 that this conclusion cannot 
be sustained for the case of gapped, frustrated chains if 
larger systems of up to 18 sites and additional values for 
the next-nearest-neighbor frustration a (a = 0.35, 0.5, 1) 
are included in the finite-size analysis. 

Additionally, several authors have recently used ana- 
lytic approaches to compute the Drude weight £>th 18 ' 20 
and the spin Drude weight D s 43 in the low-energy limit. 
If the effective model describes non-interacting particles 
such as in Ref. 18 in the case of the spin ladder, a finite 
Drude weight will naturally exist. However, to prove that 
the Drude weight is also finite in the corresponding lat- 
tice model is more difficult. The outcome of a low-energy 
description crucially depends on the effective models, i.e., 
one has to take care which operators are kept and which 
ones can be omitted when passing from the lattice to 
the continuum limit. In particular, Rosch and Andrei 38 
have shown that two independent incommensurate umk- 
lapp terms suffice to relax the spin current in massless 
models. Since this result has not been fully appreciated 
by some of the aforementioned papers, we will discuss 
the line of reasoning of Ref. 38 and apply these ideas to 
the models which are of interest in this paper. We will 
argue that a vanishing Drude weight for both kinds of 
transport is expected in generic nonintegrable (massive) 
models also. Thus, the first part of this section will be 
devoted to the discussion of transport properties in the 
continuum limit. 

In Sees. VB and VC, we will complete our numerical 
investigation of both the thermal and the spin Drude 
weight of frustrated and dimcrized spin systems with ar- 



bitrary values of a and A. The main focus will be on 
the finite-size analysis of the high-temperature prefactor 
Cth[s](-W) (see Eq. (20)). If not stated otherwise, A = 1. 

A. Bosonization 

The low-energy description of the systems studied be- 
low, i.e., dimerized chain and frustrated chains, can be 
cast in the general form of a £7(1) scalar field theory, 
known as a Luttingcr liquid (denoted by i/ix), with a 
perturbation g corresponding to a relevant operator H IC \ 
plus all irrelevant operators H m allowed by the symme- 
tries of the given problem 6 

H = i?LL + Hrel + Hin, (38) 

#ll = J dx(vK(d x ®) 2 + ^{d x 4>) 2 ) , (39) 
H Te \ = g J dxcos(acj)). (40) 

4> = 4>(x,t) is a bosonic field in 1+1 dimensions and 
is the dual field d x & = (l/K)d T <fr. K is the Luttinger 
parameter and v is the velocity. General situations in- 
volving more than one relevant operator could also occur, 
but this does not change the discussion below. 
The Hamiltonian, Eq. (38), with H lrr = corresponds 
cither to a Luttinger liquid in a generic massless situa- 
tion (g = 0, e.g., the massless regime of the frustrated 
chain) or to a sine-Gordon theory in the massive cases 
(g ^ 0, e.g., dimerized chain, massive regime of the frus- 
trated chain). These descriptions provide in general the 
correct low-energy picture if one is interested in, e.g., the 
long-distance behavior of correlators, and usually, one 
can discard the irrelevant terms H m since they only con- 
tribute with subleading corrections. 
However, as was pointed out in Ref. 38, certain operators 
have a crucial effect on transport properties and should 
therefore be taken into account to reproduce the cor- 
rect low-frequency and low-temperature behavior even 
if these operators are irrelevant in the renormalization 
group sense. The main result of Ref. 38 is that a cer- 
tain class of incommensurate umklapp operators lead to 
the decay of all currents and hence render all conduc- 
tivities finite. The emerging picture is that, except for 
very special circumstances which could happen in certain 
integrable models, one should expect a vanishing Drude 
weight and hence a finite conductivity. 
It should be stressed that Ref. 38 is devoted to massless 
cases, that is, to those situations where no relevant op- 
erators are present (g = in Eq. (40)). However, one 
can argue that the main ingredient in the proof, that is, 
the violation of all conservation laws due to the presence 
of incommensurate umklapp operators, is independent of 
the scaling dimensions of the operators involved. Hence 
one could expect a similar picture in the massive case, 
while clearly the results for k[<7](oj) will be quantitatively 
different. This conjecture is in full agreement with our 
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numerical findings. 

For pedagogical reasons let us briefly summarize the main 
results of Ref. 38: among the infinitely many irrelevant 
operators contained in H m , those which could produce a 
decay of the currents are the incommensurate umklapp 
operators, which are, in the case of pure spin models, 
generically of the form 

/ dx O nm (x) = J dxg nm cos(V2nn(j) + k nm x). (41) 

g nm are coupling constants, k nm = Inkp — mG where 
kp is the Fermi momentum, and G is a reciprocal lattice 
vector. In a fermionic representation, n is the number of 
fermions which change chirality under the action of the 
operator O n . m (x). 

These operators do not modify the low-energy expres- 
sions for the energy and spin current. The same holds 
for the relevant operator 20 in Eq. (40). The currents 
take the form (see, e.g., Rcfs. 17,65) 

-vK f 

3b = -f= / dxd x e, (42) 

V27T J 

jth = v 2 j dxd x <j>d x Q. (43) 

Comparing our notation with Ref. 38, note that j s ~ Jo 
and jth ~ Pt- The spin current neither commutes with 
the relevant operator nor with H lvr while for the thermal 
current, [H leh jth] = but [H ilr , j th ] ^ 0. 
The key observation 38 is that in the presence of one such 
operator C n m (a;), there is still a conserved current which 
can be written as a linear combination of the spin and 
thermal current 

^conserved k nm j s + 2nj th . (44) 

However, as soon as more than one of the operators 
On,m{ x ) a re considered, no conservation law of the type 
(44) survives and hence the conductivity is expected to 
be finite. Since there is no reason a priori to exclude such 
incommensurate operators, this seems to be the generic 
situation. As incommensurate operators have been con- 
sidered neither in Ref. 18 nor in Ref. 20 it is not clear 
whether their results of a finite thermal Drudc weight in 
the low-energy limit provide a proof of D t h{T > 0) > 
for the respective nonintegrable lattice models. 

Generally, even if the current operator is not conserved, 
a nonzero Drudc weight can be caused by the existence of 
(nontrivial) conserved quantities {Qi} with a finite pro- 
jection on the current operator jth[s] m the Liouvillc space 
(sec, e.g., Rcfs. 38,46,47,57). More precisely, the Drude 
weight is nonzero if 

Ai[s] (T) = jsJ^ttthM I Pjth[.]) > (45) 

where V is the projection operator on all conserved quan- 
tities {Qi}- (A | B) denotes Mori's scalar product 57 in the 



space of operators 

(A(t)\B) = ± J dr(A(^B(ir)). (46) 
o 

Under certain circumstances, e.g., integrability, it is pos- 
sible to construct an infinite set of {Qi} (see Ref. 30 and 
references therein). Still, the evaluation of Eq. (45) is a 
nontrivial problem. 

In the literature 30,43 , one often refers to a weaker condi- 
tion than Eq. (45), namely Mazur's inequality 30,66 where 
a subset of the {Qi} or even only one operator Qi £ {Qi} 
is taken into account. Therefore, Mazur's inequality pro- 
vides a lower bound for the Drude weight. For instance, 
the conservation law 38 , Eq. (44), would suffice to prove 
that D s > if no further incommensurate operators are 
considered. 

Quite recently, such an operator Qi has been found for 
charge transport in a Luttingcr-liquid plus interactions 
spoiling both the integrability and the conservation of 
the spin current operator 43 , reformulating Zotos and co- 
workers^ results from Ref. 30 in the continuum limit. 
The conserved quantity can be readily identified as the 
thermal current operator (compare Refs. 17,20). How- 
ever, their proof of (j s \Qi) > assumes particle-hole 
symmetry to be broken, e.g., by the existence of a mag- 
netic field (see also Refs. 22,30). Consequently, for the 
class of models considered in Ref. 43, a nonzero Drude 
weight can be inferred. Since the incommensurate oper- 
ators of Eq. (41) are explicitly excluded in Ref. 43, their 
result does not contradict our numerical indications for a 
vanishing Drudc weight D s in nonintegrable spin-lattice 
models. 



B. Thermal transport in nonintegrable models 

Frustrated chain - In the thermodynamic limit, 
the low-energy spectrum of a frustrated chain with 
a < a cr it ~ 0.241 67 is gaplcss and gapped for a > a cr it- 
The thermal Drudc weight in the gapped regime of 
frustrated chains has been discussed in detail in Rcfs. 
17,21 where we found clear indications of a vanishing 
Drude weight for N — ► oo. Figure 8 shows the thermal 
Drude weight D t h(N, T) and the specific heat Cy for 
a = 0.2 and N = 8, 10, 12, 14, 16, 18. For chains of finite 
length, the data at low temperatures are dominated by 
the finite-size gap. Hence the Drude weight and specific 
heat are exponentially suppressed for small T. While 
the specific heat converges to the thermodynamic limit 
at temperatures T > 0.25J, strong finite-size effects 
are present in the data for the Drude weight at all 
temperatures. 

At low temperatures, D t h{N,T) monotonically in- 
creases with system size similar to the case of a = 0.35 
(see Fig. 3 in Ref. 17). In Ref. 21, we have argued 
that the notion of an increasing Drude weight at low 
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A=l, a=0.2 (N=8,...,18) 

— ED: ° lh 
ED: C v 

□ HTSE: C„ 




FIG. 8: Thermal transport, frustrated chain: Drude weight 
D t h(N,T) (dashed lines) and specific heat Cv (solid lines) for 
a = 0.2 and TV = 8, 18 sites (arrows indicate increasing 
system size). The plot includes data for Cv from a high- 
temperature series expansion (HTSE) reproduced from Ref. 
68 (note that only the bare high-temperature series up to or- 
der 10 in J/T is shown. By means of extrapolation schemes, 
the HTSE can be extended to significantly lower tempera- 
tures, see Ref. 68). 



temperatures does not support the conjecture 15 of a 
finite D t h for TV — > oo for a = 0.35. In fact, we have 
shown that a crossover temperature T* which wc define 
by D th (N + 2,T*) = D th (N,T*) and even TV seems to 
extrapolate to zero as a function of system size. This 
implies that the temperature range where one observes 
an increasing Drude weight with system size could 
vanish for TV — > oo. An analogous finite-size analysis 
of T* for a = 0.2 (not shown in the figures) could also 
be interpreted in the same sense, i.e., the temperature 
interval where D th (N,T) < D th (N + 2,T) tends to 
vanish for TV — > oo. To summarize the discussion of 
the low-temperature regime, we emphasize that the 
finite-size data should not be used to speculate about 
the thermodynamic limit. 

In Fig. 9 (a), we present the high-temperature 
prefactor Cth (AO for several values of a both in the 
gapless (a = 0.1,0.2, open symbols) and the gapped 
regime (a = 0.25,0.35,0.5,1, solid symbols; the last 
three sets have already been shown in Ref. 21). Finite 
systems with up to TV = 18 sites have been analyzed. 
While Cth (TV) appears to be almost constant in the case 
of a = 0.1, a substantial decrease with system size is 
observed for larger a and sufficiently large TV as it is 
especially obvious for a = 1. We also note that, regard- 
ing thermal transport, the data at the Majumdar-Ghosh 
point a = 0.5 (see Ref. 69 and references therein) do not 
point to any peculiarities. 

In panel (b) of Fig. 9, C t h(TV) is plotted versus a for 
TV = 14, 16, 18. Starting at small a, we observe that 
C t h (TV) is discontinuous at a — which will be com- 
mented below. The curve further decreases with a and 
exhibits a minimum atas 0.4 for TV = 14 and a ss 0.5 
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FIG. 9: Thermal transport: high-temperature prefactor 
Cth (TV) for frustrated chains. Panel (a): Cth (TV) /Cth (8) ver- 
sus 1/TV for various values of a both in the gapless and gapped 
regime (TV = 8, 18) (see also Ref. 21). Panel (b): C th (TV) 
for TV = 14, 16, 18 as a function of frustration a. The arrow 
indicates the system size independent value of Cth at a = 0. 



for TV = 16. The position of the minimum seems to be 
further shifted towards larger a on growing TV. Further 
increasing the next-nearest-neighbor interaction drives 
the system into the limit of two decoupled chains each 
with TV/2 sites and interchain interaction a. Exactly 
for J — > 0; Ja = const, the current operator is again 
conserved. Consequently, one expects the Drude weight 
to increase for large a at finite and fixed TV. This feature 
is indeed found for a > 0.5, see Fig. 9 (b). 

Fig. 9 (b) indicates a difference between the gapped 
and gapless regimes: the decrease of C t h{N) with TV is 
weaker in the gapless regime. We suggest the following 
two scenarios for further discussion: (i) the Drude 
weight is nonzero in the gapless regime and zero in 
the gapped regime; (ii) the Drude weight is zero for all 
a > 0, but, depending on a, there is a characteristic 
system size N(a) with Cth(TV) w const for TV < N(a) 
and monotonically decreasing for TV > TV (a). 

The first interpretation might be plausible in view of 
the significant differences in the low-energy properties 
for a < a cr it and a > a CT it- However, since Cth is 
essentially the Drude weight at infinite temperature 
where all states contribute with equal weight, it is not 
clear why low-energy features should play a crucial role 
for the finite-size scaling in the limit f3 — » 0. 
A second objection against the first scenario arises from 
the analysis of the level-spacing distribution both in the 
gapped and gapless regime. Exploiting translational 
invariance and conservation of total S* ot already lifts 
all degeneracies on finite systems as is obvious from 
Fig. 10 showing the integrated level spacing distribution 
/(e); sec Eq. (32). The difference to the spectrum of 
the integrable model (see Fig. 3) is striking: while a 
large fraction of states with AE n < 10~ 8 J is present 
for A = 0.5, a = 0, no such candidates for degenerate 
states appear in the case of cv = 0.2, A = 1. This 
feature is characteristic for a > which, in particular, 
supports the conjecture that transport properties in the 
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FIG. 10: Frustrated chain: distribution of level spacings in 
the spectrum of finite chains with 8 < iV < 18 (bottom to top 
as indicated by the arrow) for a = 0.2. AE„ is the difference 
of adjacent energy levels in subspaces classified by total Stot 
and momentum k. The number 7(e) of AE„ with AE n < e 
summed over all subspaces is plotted versus e. 



gapped and gapless regimes should not be different at 
high temperatures. Exceptions are found for a = 0.5, 1. 
At the Majumdar- Ghosh point, one degenerate state 
occurs if N/2 is even. In the latter case (i.e., a = 1), 
there are degenerate states in the spectra of chains with 
N = 10, 12, 14, 18 which are, however, small in number 
(w 10 for N = 18). 

Next, there is the discontinuity of C t h(N) at a = 0. 
A small, but finite frustration (e.g., a = 10~ 4 , 10~ 3 ) has 
the effect that degeneracies are lifted while the values of 
diagonal matrix elements |(n|j s |n)| 2 are almost unaf- 
fected. This leads to the substantial difference between 
Cth at a — compared to small, but finite a > 0. 
Finally, we mention that the fact of Cth(N) const for 
small a could be a consequence of the proximity to the 
integrable point a = 0. 

In conclusion, our numerical data indicate a vanishing 
thermal Drude weight for arbitrary values of a at high 
temperatures. This result is difficult to reconcile with 
the recent findings 18,20 of a nonzero Drude weight in 
the continuum limit as we have discussed above that a 
crossover from a nonzero to a zero Drude weight as a 
function of temperature is not likely. 



Spin ladder and dimerized chain - Now we turn to the 
cases of the dimerized chain and the spin ladder. First, 
we discuss the numerical data for the thermal Drude 
weight D t h(N,T) taking the example of a spin ladder 
with a = Jy/Ji = 0.5, J± = J for N = 8, 10, 12, 14, 16. 
Second, the results from a finite-size analysis of the high- 
temperature prefactor for both spin ladders and dimer- 
ized chains are presented. Finally, we comment on pos- 
sible implications of our results for the interpretation of 
recent experiments 3 . 

Due to the dimerization, the unit cell of our model is 
doubled restricting the maximal system size to N = 16 



0.5 




FIG. 11: Thermal transport, spin ladder (a = 0.5, A = 0, 
N = 8, ... , 16). Panel (a): specific heat CV and susceptibility 
X- Panel (b): Drude weight D th (N,T). ED for D t h,C v ,X- 
lines; QMC for CV(x): squares(circles) reproduced from Ref. 
15(70). 



in our numerical computations at present. 

Regarding the level-spacing distribution, degeneracies 
are still present in the case of the spin ladder. For ex- 
ample, there are < 10 2 degenerate states for N = 16 
compared to > 10 3 in the integrable case (see Fig. 3). 
The spectra of dimerized chains show the same features. 

The specific heat Cy, the susceptibility x, and the ther- 
mal Drude weight D t h(N, T) are plotted versus T in Fig. 
11 for N = 8, 10, 12, 14, 16 and a = 0.5, A = (panel (a): 
Cy, x, panel (b): D t h(N,T)). The main characteristics 
are the following: (i) for the specific heat, finite-size ef- 
fects are small and negligible for the susceptibility; (ii) 
the data for Dt^NjT) display strong finite-size effects 
at all temperatures; (iii) Dth(N,T) is monotonically de- 
creasing at high temperatures for N > 8 and T > 0.6J; 
(iv) the positions of the maxima of the specific heat and 
the Drude weight are different; (v) for N/2 even(odd), 
the data are monotonically increasing(decreasing) at low 
temperatures. The latter may be attributed to the fact 
that D t h(N, T) is diverging for decoupled chains and odd 
N. 

It should be stressed that the restricted number of sys- 
tem sizes analyzed here precludes any conclusions from 
the finite-size scaling at temperatures T < 0.6J; in par- 
ticular, since the iV-dependence is nonmonotonic. Note 
that even a monotonic increase of D t h (N, T) with system 
size at low T as observed in the case of frustrated chains 
with, e.g., a = 0.35 does not unambigously point to a 
finite Drude weight (see the discussion of D t h(N,T) of 
frustrated chains and Ref. 21). 

While our numerical data for the specific heat are in 
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FIG. 12: Thermal transport: Cth(-W) f° r spin ladders (panel 
(a)) and dimerized chains (panel (b)). Data are shown for 
TV = 8, . . . , 16 spins (i.e., ladders with 4, . . . , 8 rungs). 



quantitative agreement with the results from Ref. 15 for 
the same choice of parameters, only qualitative consis- 
tency is found for the thermal Drude weight regarding 
points (ii)-(iv). Note that the Drude weight Dth in Ref. 
15 is measured in units of Jii ~ a J instead of J± = J 
used in this work. The specific heat, however, is di- 
mcnsionlcss. Also, compared to point (v), the opposite 
monotony behavior is observed at low temperatures in 
Ref. 15. One may speculate that this difference is due 
to the use of twisted boundary conditions and different 
definitions for the energy-current operator j t h in Ref. 15. 

The finite-size analysis of the high-temperature pref- 
actor Cth(-W) reveals a systematic decrease with system 
size for N > 8 both for the case of the spin ladder and 
dimerized chain as it is evident from Fig. 12. In particu- 
lar, by normalizing the data on the respective values for 
Cth(N = 8) the finite-size dependence of the data for the 
spin ladder (panel (a) in Fig. 12) appears to be almost 
independent of the interchain coupling a = J\\/ J± for the 
choice of parameters considered here (a = 0.1,0.5,1,2) 
including the isotropic ladder (a = 1). Only the results 
for the dimerized chain with A = 0.9 show less evidence 
for a vanishing of Cth- This docs, however, not ques- 
tion the conclusion of a vanishing thermal Drude weight 
because A = 0.9 is still very close to the homogeneous 
Heisenberg chain where D t h(T > 0) is finite. 

To summarize the finite-size analysis, one can con- 
clude that the numerical data for D t h of spin ladders and 
dimerized chains indicate a vanishing Drude weight for 
N — > co. In particular, this includes the isotropic spin 
ladder (Jj_ = Ju;a = 1) which is of relevance because 
the magnetic properties of LasCa9Cu2404i are well de- 
scribed by J_l w J|| 71,72 . 

Recently, first attempts have been made to extract 
magnetic mean free paths Z mag from the experimental 
data for the magnetic part re mag of the thermal conduc- 
tivity of La5CagCu2404i. Assuming that heat is car- 
ried mainly by the elementary excitations of spin ladders 
(i.e., dispersive triplet modes) Hess et al. 5 have used a re- 
laxation time ansatz for the respective kinetic equation 



k 



Cv,fc Vk Ik- 



(47) 



Here, CV,fc is the specific heat per fc-space volume, Vk is 
the triplet dispersion, and Ik the momentum-dependent 
mean-free path. This approach results in very large 
mean-free paths Z mag of the order of 3000A at 100K cor- 
responding to w 770 lattice constants. Following this 
work, Alvarez and Gros 15 have suggested a much smaller 
value for Z mag , namely, 176 A as 45 lattice constants at 
T = 100/\ . They have attributed this significant differ- 
ence to the fact that D t h(N,T) <fi Cy as is seen in the 
numerical data. It is, however, straightforward to check 
that Eq. (47) leads to K mag ^ Cy, implying that this 
alone does not explain the different values found for £ mag 
in Refs. 3,15. 

The result for ^ mag of Ref. 15 is explicitly based on the 
assumption of a finite Drude weight D t h for spin ladders 
which is questionable in view of the detailed numerical 
results presented in this paper. 



C. Spin transport in nonintegrable models 

In this last section, we give an overview of our results 
for the Drude weight for spin transport in nonintegrable 
models. The dependence of D s on frustration and 
dimerization will be systematically discussed. Our 
analysis of the Drude weight at finite temperatures 
includes next-nearest-neighbor interactions J2i Si ■ S1+2 
extending previous numerical studies of nonintegrable 
lattice models 29,33,41 where different kinds of Ising-likc 



interactions Sf+v * = 2,3) have been consid- 

ered. Note in this context that frustration cannot be 
treated with QMC simulations due to the sign problem 41 . 

Frustrated chain - Results for the Drude weight 
D^' II (N, T) are shown in Fig. 13 for a = 0.2 (main panel) 
and a = 0.5 (inset). In the gapless regime (a < a cr it), 
the finite-size data display features similar to the XX Z 
model: (i) D^(N,T = 0) > 0; (ii) D* T (N, T) w const 
at small temperatures; (iii) D*(N,T) ~ D^ T (N, T) at 
high temperatures, but significant deviations at low 
temperatures; (iv) a monotonic decrease with system 
size at high temperatures T > 0.5J. Similar to the 
case of the thermal Drude weight of frustrated chains, 
we observe that for even N, the monotonic decrease 
of D^' 11 (N,T) with system size at high temperatures 
changes to an increase at lower T. However, the tem- 
perature, where this change in the monotony behavior 
occurs, is strongly shifted to lower temperatures as N 
grows; see Fig. 13. This resembles the case of Dth{N,T) 
for a = 0.35 discussed in Ref. 21 and we conclude that 
the numerical data for D s at low temperatures do not 
give unambiguous evidence for a finite Drude weight. 

The observation of 1 (N, T = 0) > for small a has 
also been reported in Ref. 26. There, using the Lanczos 
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FIG. 13: Spin transport, frustrated chain: Drude weight 
Di'"(N : T) for N = 8, 11, 12, 13, 15, 16, 18 and a = 0.2, A = 1 
(dashed lines: D a 7 (iV,T); solid lines: D^(N,T)). Arrows in- 
dicate increasing systems size. Inset: Drude weight at the 
Majumdar-Ghosh point (a = 0.5, A = 1) for N = 12, 16. 
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FIG. 14: Spin transport: high-temperature prefactor C S (N) 
for frustrated chains. Panel (a): C S (N)/C S (8) versus 1/N for 
various values of a both in the gapless and gapped regime 
(N = 8,9, 18). Panel (b): C S (N) for N = 14, 16, 18 as a 
function of frustration a. The inset shows a blowup for small 
a and N = 14, 16. 



method and truncation in the S^ ot = subspace, a 
nonzero Drude weight at T — has been found for 
a < 0.43 and N = 20. To clarify whether D B (T = 0) > 
survives in the thermodynamic limit, one should exploit 
Kohn's formula for T = using, e.g., the Lanczos 
algorithm, which is, however, not the purpose of the 
present paper. 

Consistent with Ref. 26, we find D^(N,T = 0) = 
in the case of a = 0.5 (see inset of Fig. 13). Notice that 
D*(N,T) ~ D**(N,T) at all temperatures for N = 16. 
We now turn to the question of a nonzero Drude weight 
in the thermodynamic limit by a finite-size analysis of 
the high-temperature prefactor C S (N). C S (N) is plotted 
versus 1/N in Fig. 14 (a) for a = 0.2,0.35,0.5,1. First, 
C B (N) monotonically decreases with system size (except 
for odd-even effects) for all values of a presented here 
and exhibits a discontinuity at a = analogous to 
Cth{N) (see inset of Fig. 14). Second, the data for 
a = 0.2 may in principle be extrapolated to a finite 
value in the thermodynamic limit (see, however, remarks 
below). The behavior for a > a cr it is similar to the 
case of the thermal Drude weight since C S (N) decreases 
rapidly with N and faster than 1/N. 

Another remarkable difference between the gapped 
and the gapless regime is revealed in Fig. 14 (b) where 
we show C S (N) as a function of frustration a for 
A" = 14, 16, 18. In contrast to the thermal Drude weight 
(see Fig. 9), C S (N) first grows with a and exhibits a 
maximum around a ~ 0.2. For larger a > 0.7, C S (N) 
increases with a again analogous to the case of the 
thermal Drude weight (see Fig. 9). 

We interpret the data for a > a cr ;t in terms of a 
vanishing Drude weight for A^ — > oo. For the gapless 
regime, the possibility for a nonzero C s cannot be ruled 
out by our data although the absence of degeneracies 
(see Fig. 10) supports the conclusion of C s — for all 
a > 0. 
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FIG. 15: Spin transport, dimerized chain: Drude weight 
Dl' n (N,T) for TV = 10, 16 (top to bottom as indicated 
by the arrow) and a = 0,A = 0.5 (dashed lines: D' S (N,T); 
solid lines: D^ 1 (N,T)). In the inset, the exact result from 
Eq. (50) (dot-dashed line) for the Drude weight of a dimer- 
ized XY chain with A = 0, a = 0, A = 0.5 and numerical data 
(solid line) for N — 16 sites are shown. 



Dimerized chain, spin ladder - While for the frustrated 
chain and the spin ladder the fermionized Hamiltonian, 
Eq. (5), contains interaction terms even at A = 0, the 
case of the dimerized XY model, (A = 0, A ^ l,a = 
0) corresponds to a model of free, but massive fermions 
which can be solved exactly (see, e.g., Ref. 18). We will 
start with a discussion of this limiting case where a finite 
Drude weight D s exists. 

The Hamiltonian in terms of spinless fermions (see Eq. 
(5)) reads 

^ y = Ey( c U+H.c) (48) 
l 

(A; = A for I even and A/ = 1 otherwise). A straightfor- 



17 



(a) Dimerized chain 
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FIG. 16: Spin transport: high-temperature prefactor C S (N) 
of the Drude weight for the dimerized chain (panel (a)) and 
the spin ladder (panel (b)); both for N < 16. 



ward computation diagonalizes H XY (see, e.g., Ref. 18 
for details) 

H XY = J> (<+%+ - 4,-°*,-)' (49) 



leading to two modes with a gapped dispersion = 
J \J (1 — A) 2 /4 + Acos(fc). Obviously, the spin-current 

operator j s — J2k v k( a l,+ a k,+ ~~ a l.- a k,-) ' 1S conserved 
and the Drude weight can be computed exactly 



D -< T > - w I 



dk- 



cosh z [ efe /(2T)] 



(50) 



Vk = dek/dk is the velocity. 

In Fig. 15, the Drude weight on finite systems is 
plotted for a dimerized chain with A = 0.5 and both 
A = (inset) and A = 1 (main panel). In the former 
case, we see that the numerical data for 16 sites agree 
with the exact expression for N — > oo, Eq. (50), at 
temperatures T > 0.2J. The same holds for smaller 
system sizes (not shown in the figure) for high T. This is 
in contrast to the curves for A = 1 (main panel), where 
no convergence is observed at all temperatures. 
For A = 1, both D*(N,T) (dashed lines) and D* T (N,T) 
(solid lines) have been evaluated proving their equiva- 
lence at high temperatures. Small systems (N = 10, 12) 
still exhibit a large nonzero value of Df 7 (7V, T = 0) 
which rapidly decreases with system size. 

The data for the high-temperature prefactor C S (N) 
are collected in Figs. 16 (a) (dimerized chain) and (b) 
(spin ladder). A substantial decrease of C S (N) with 
system size is observed for all choices of parameters and 
A > 8, indicating a vanishing Drude weight at high 
temperatures. 

In summary, our numerical data yield evidence for 
nonballistic spin transport in nonintegrable spin models 
with a spin gap, i.e., D S (T > 0) = in the thermody- 
namic limit. This result confirms Zotos and coworker's 



original conjecture 29 that transport in nonintegrable 
models should be normal. The system sizes investigated 
here may, however, be too small to clarify whether the 
Drude weight is zero or not in the gapless phase of the 
frustrated chain. 



VI. CONCLUSION 

In this paper, we have presented a detailed study of 
the Drude weights for thermal and spin transport in one- 
dimensional spin systems at finite temperatures. Let us 
focus here on the main results while more detailed sum- 
maries can be found in the preceding sections. 

Thermal transport in the XX Z model is generically 
dissipationless due to the conservation of the energy cur- 
rent operator 30 . Our data for the thermal Drude weight 
are in excellent agreement with the Bethe ansatz 14,49 . 
For spin transport in this model, we have presented nu- 
merical results for various values of the anisotropy. Our 
data confirm the observation 29,32 ' 33, 37,41,43,44 £ a g n j^ e 
Drude weight D s in the gapless regime (|A| < 1). We 
have discussed some so far unresolved issues: first, the 
exact temperature dependence of the Drude weight in 
the critical regime (|A| < 1) and, second, the question 
of whether D S (T > 0) is finite for the isotropic chain 
(A = 1). Regarding the first point, analytical 43 and nu- 
merical results 41 are compatible with the finite-size data 
presented here. Regarding the second issue, the exact 
diagonalization data of finite systems with N < 18 do 
favor a finite Drude weight at A = 1 . 

In the case of nonintegrable models (frustrated chain, 
dimerized chain, spin ladder) , our main result is that the 
finite-size analysis of the ED data does not indicate a 
finite Drude weight in the thermodynamic limit either 
for thermal or spin transport, but rather supports the 
conclusion that transport in these systems is dissipativc. 
While we have concentrated our numerical analysis on 
the finite-size scaling at high temperatures, this result 
is corroborated by bosonization in the low-temperature 
limit. 
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